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The problem of the quantum-limited or intrinsic linewidth 
of a good-cavity laser is revisited. Starting from the Scully- 
Lamb master equation, we present a fully analytical treatment 
to determine the correlation function and the spectrum of the 
cavity field at steady state. For this purpose, we develop an 
analytical approximation method that implicitly incorporates 
the microscopic fluctuations of both the phase and intensity of 
the field, and, in addition, takes full account of the saturation 
of the nonlinear gain. Our main result is a simple formula for 
the quantum-limited linewidth that is valid from near to far 
above threshold and also includes the presence of thermal pho- 
tons. Close to the threshold, the linewidth is twice as large 
as predicted by the standard phase-diffusion treatment ne- 
glecting intensity fluctuations, and even 50% above threshold 
the increase is still considerable. In general, quantum fluctu- 
ations of the intensity are present and continue to influence 
the linewidth as long as the photon-number distribution is not 
strictly Poissonian. This inherent relationship is displayed by 
a formula relating the linewidth and the Mandel Q-parameter. 
More than 100% above treshold the linewidth is found to be 
smaller than predicted by the standard treatment, since the 
simple phase-diffusion model increasingly overestimates the 
rate of phase fluctuations by neglecting gain saturation. In 
the limit of a very large mean photon number the expected 
perfectly coherent classical field is obtained. 
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I. INTRODUCTION 

When the radiation field is described classically and 
spontaneous emission is neglected in comparison to in- 
duced emission, the steady-state linewidth of an ideal 
single-mode laser, i. e. of a laser with a perfectly sta- 
bilized classical intensity, has the shape of a J-function. 
This is due to the fact that the cavity losses are exactly 
compensated for by the gain and, as a result, the field in 
the resonator remains constant. From a fully quantized 
description of the electromagnetic field, however, it fol- 
lows that the laser linewidth cannot be smaller than a cer- 
tain quantum limit, related to the well-known Schawlow- 
Townes linewidth, which is inversely proportional to the 
laser intensity [Q. This limit has been studied exten- 
sively in the past decades (see Refs. and references 
therein) . In view of the importance of stable coherent sig- 
nals for various high-precision measurements and because 
of an ongoing interest in fundamental problems of quan- 



tum optics, renewed attention has been paid recently to 
the quantum limitation of the laser linewidth. The inves- 
tigations have been extended to cover bad-cavity lasers 
0, lasers without inversion and a chaotic laser cavity 
[p|. Different systems have been devised in this context 
to reduce the quantum limit of the linewidth by means 
of correlated spontaneous emission schemes [|| . Recently 
even theoretical models for light amplification without 
stimulated emission have been investigated in order to 
obtain a reduced ultimate quantum limit of the laser 
linewidth (lO). 

The quantum-limited laser linewidth is also called the 
intrinsic or natural linewidth of the laser. Its origin, like 
the origin of the microscopic intensity fluctuations, lies 
in the fact that in the stationary regime of operation the 
balance between the gain and loss processes sustains a 
constant average field only but, due to the quantum na- 
ture of these processes, fluctuations of the field around its 
mean are induced on a microscopic scale. We note that 
the resonator losses are caused by the outcoupling of the 
field through the output mirror as well as by any addi- 
tional linear damping process such as absorption. For 
good-cavity lasers the combined effect of these losses can 
be described by a single cavity-damping constant 7 that 
is the sum of the constants referring to the individual 
processes. The usual treatment of the intrinsic laser 
linewidth rests on the approximation that the linewidth 
arises from fiuctuations of the phase of the field described 
by phase diffusion. 

In this paper we rely on the Scully-Lamb model |^ 
of the laser since it is applicable for an arbitrary strong 
saturation of the lasing atomic transition. Neglecting in- 
tensity fiuctuations, a simple analytic expression for the 
linewidth has been derived in this model by means of 
different methods |^-||. On the other hand, the intrinsic 
laser linewidth can be determined exactly by numerically 
calculating the first-order correlation function of the field 
and performing the Fourier transform of the latter. Thus 
both the effects of fluctuations of the phase of the fleld 
and of its amplitude, or intensity, respectively, are im- 
plicitly taken into account. Investigations of this kind 
have been performed by N. Lu JO who started from the 
Scully-Lamb master equation ||^and found numerically 
that near threshold the intrinsic laser linewidth is up to 
twice as large as that given by the phase-diffusion coeffi- 
cient at the mean photon number. 

In the present contribution we derive an analytical 
expression for the quantum-limited laser linewidth by 
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means of investigating the two-time first-order correla- 
tions of the field. The treatment is restricted to single- 
mode lasers in the good-cavity limit, i. e. we make the 
usual assumption that the cavity damping time 7"^ is 
long in comparison to all other relevant time scales. This 
ensures that the time dependence of the polarization and 
population inversion of the gain medium can be adiabat- 
ically eliminated and the Scully-Lamb model of the laser 
can be applied. Our approximation scheme is an exten- 
sion of an analytical method developed previously by one 
of the authors for calculating photon-number variances 
| p2| , p^ . It makes use of the fact that the photon- number 
distribution of the laser radiation is strongly peaked at 
a large mean photon number. Therefore it is not nec- 
essary to study the density-matrix elements of the field 
in detail but it suffices to directly evaluate the expecta- 
tion values and correlation functions of interest, in the 
approximation of small fluctuations. 

The paper is organized as follows. In Sec. II we de- 
velop a general approximation method for the determi- 
nation of the first-order correlation function of the field 
in the resonator of a micromaser or laser. The method 
is applied in Sec. Ill to study the laser linewidth. In 
order to reveal the influence of intensity fluctuations on 
the latter, an expression for the photon number variance 
of the laser in the presence of thermal photons is also 
derived in this section. Our results are discussed in Sec. 
IV and compared to the standard linewidth formula en- 
suing from the phase diffusion approximation. Finally a 
summary is given in Sec. V. 



II. GENERAL APPROXIMATION METHOD 

The power spectrum of a radiation field is defined as 
the Fourier transform of its normalized first-order cor- 
relation function. When the field is a single-mode field 
with frequency the steady-state spectrum is given by 



S[uj) = -Re / 



(2.1) 



where a and are the photon annihilation and creation 
operators in the interaction picture and the subscript s 
denotes the stationary state. We consider a single-mode 
radiation field contained in a leaky cavity and being sus- 
tained by a gain mechanism. The overall losses of the 
field mode are supposed to be due to the coupling of the 
cavity field to the environment, modeled by a reservoir in 
thermal equilibrium. Therefore the damping of the field 
can be described by the standard master equation p = Lp 
for its reduced density operator p, where the action of the 
superoperator L is defined as 



7 ' ' + 

Lp = "■^(1 + ni,){a^ap — 2apa^ + pa' a) 

— ^nb{aa^p — 2a) pa + paa^). 



Here 7 is the cavity damping constant, and denotes 
the mean number of thermal photons in the cavity. We 
also note that for any field operator a the relations. 



Tr(aticr) = -iTr(aV), 
Tr(ataLo') = — 7Tr(a^a(T) -f 7?iftTrcr, 



(2.3) 
(2.4) 



follow from Eq. (2.2). To model the gain in a simple way, 
we assume that excited atoms are injected into the cav- 
ity with rate r. The atoms are supposed to interact with 
the field one after the other during a time r that is neg- 
ligibly short in comparison to both the cavity damping 
time 7^^ and the mean time interval between successive 
atoms. The effect of a single atom on the field can be 
formally written as p[t + t) — Mp{t), where the super- 
operator M depends on the specific kind of interaction 
process. For a micromaser the interaction time r is given 
by the transit time of the atoms through the microwave 
cavity, and M = M{t) has to be obtained in the usual 
way from the Jaynes-Cummings Hamiltonian p5| for 
the atom-field interaction by performing the trace with 
respect to the atoms. In order to describe a laser, it is 
assumed that excited two-level atoms are injected into 
a resonant cavity and interact independently with the 
field during time intervals that are determined by the 
survival time F"^ of the atoms as effective two-level sys- 
tems, before they decay into different energy states [Q. 
The superoperator M that has to be used for the laser is 
found by averaging M(r) with respect to the interaction 
time r, according to M = F M(r)e-^^dT ||]. 

When the injection times of the atoms are uncorre- 
lated, i. e. for Poissonian pumping, the evolution of 
the field due to the combined action of the gain and loss 
mechanisms obeys the master equation 

p^r{M -l)p + Lp, (2.5) 

which has the formal solution p{t) ~ V{t)p{0) with 

y(t) = cxp{[r(M- 1) + L]t}. (2.6) 

Due to the Markovian character of the master equation 
the two-time correlation function necessary to determine 
the spectrum can be easily calculated. For the stationary 
state we find 



(at(t)a(0)). = Tr[at^(i)(ap)], 



(2.7) 



(2.2) 



with p = limt^oo p{t) denoting the steady-state density 
operator. Making use of Eq. (2^) and taking int o ac- 
count the relation ( |2.3[ ), we obtain from Eq. (2.7) the 
differential equation 

^(at(t)a(O)), = [r[b{t) - 1] - ^} (at(t)a(0))„ (2.8) 
where we have introduced the abbreviation 
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bit) 



TT[a^M<j{t)] ^1 



(at(t)a(O)), 



with 



(2.9) 



(2.10) 



Under steady-state conditions it is possible to eliminate 
the injection rate r from Eq. (2^) by expressing it in 



terms of field expectation values and the cavity decay 
rate. To this end we start from the steady-state relation 



As required, the linewidth reduces to the empty-cavity 
linewidth 7 when r = or when b = 1, i. e. when 
either no atoms are injected at all, or when the atoms 
do not interact with the field and M is equal to the unit 
operator. Under these conditions the steady-state is, of 
course, the thermal state with mean photon number nt, 
or, for rifc — 0, the vacuum. 

We assume that the superoperator M in Eq. (2.5) is 



due to resonant one-photon interaction of the field with 
a two-level atom, initially in its excited state. In the 
photon-number representation, M has the general form 



(2.16) 



^n = Y^ np^ .^ = Tr[a}aVp\ = 0. 

n=0 



(2.11) 



By inserting Eq. (2^) and utilizing Eq.(2^) we arrive at 
the photon-number balance equation 



[YT{a) aM-p) 



-l{n-nb), 



(2.12) 



where n = {a^a)s = Tr(a^ap) is the steady-state mean 
photon number (a^(O)a(O)). 

So far all equations hold exactly and can be applied to 
lasers as well as to micromasers with Poissonian pump- 
ing. In order to obtain an analytical expression for the 
spectr um , it is necessary to calculate the quantity b in 
Eq. ( ^.8[ ). For this purpose we use an approximation 
method that rests on the assumption that the steady- 
state photon-number distribution in the cavity is strongly 
peaked at a large mean photon number n. Since the or- 
der of magnitude of the width of the photon-number dis- 
tribution is determined by vAr?, we assume that the 
relations 



n > 1, 



< n 



(2.13) 



are fulfilled for the mean photon number and its variance 
An^ . By means of performing an expansion with respect 



to suitable parameters in Eq. (2.9) and replacing a{t) 
by its initial value tT(0) in small terms in this expansion, 
it is possible to obtain an approximate expression for b. 
It will not depend on time, provided that the leading 
term in the expansion proves to be time-independent. In 
this case the value of b depends only on p and on the 
specific form of M, i. e. the specific kind of the atom- 
field interaction process, and we can easily integrate Eq. 
(O) to obtain 



(at(t)a(O)). 



,[r(fc-l)-^]t 



(2.14) 



Fourier transformation according to Eq. (2.1) yields a 
Lorentzian steady-state spectrum that is centered at the 
frequency ujf) = v ^Im6. The linewidth (full width at 
half maximum) is given by 



Atj = 7 + 2r(l-Re6). 



(2.15) 



the cases of a micromaser or a laser, respectively. It then 
follows that the steady-state density operator is diagonal 
in the photon-number representation, as ca n be easily 
shown with the help of Eqs. (2^) and (|2.5[ ). Therefore 
the expectation value of the field (which in the steady 
state corresponds to its time average) vanishes. 



Tr(ap) = 0. 



(2.17) 



In terms of the amplitude an d ph ase of the expectation 
value of the field, (a)s, Eq. ( ^.17j) requires that for pa- 
rameter values for which the amplitude is fixed, the phase 
is uniformly distributed between and 27r, i. e. all phase 
values are equally probable at steady state. We note that 
this holds for lasers as well as for micromasers unless the 
atoms are injected in a definite superposition of their en- 
ergy states. 

In general, a correlation exists between the values of 
the field at different times that decays with increasing 
time difference. It is this decay that deter mines the 
spect rum and the linewidth according to Eqs. (2.1) and 

(iH. 



III. THE LASER LINEWIDTH 



For the case of the laser, the coefficients An^m and i?„,m 



in Eq. ( 2.16 ) can be written as (|-||] 

_ X [1 + \{n + m)] + \x^{n - mf 
2 + x{n + m + 2) + \x{n-mY 

and 



2 -|- x(n -f m -|- 2) + ■aX^in ~ m) 



(3.1) 



(3.2) 



where we introduced the saturation parameter x = 
4g^/r^ with g and F"^ denoting the atom-field coupling 
constant and the average lifetime of the atom as a two- 
level system, respectively. After cha ngin g the index of 
summation appropriately, from Eqs. ( ^.9[ ) and ( ^.16 ) we 
obtain 
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b{t) = 1 + 



^ n— 1 



\/ncr„_i,„(i) 



(3.3) 



Here use has been made of the fa ct th at the terms pro- 
portional to in Eqs. (3J^) and ( |3.2| ) can be neglected 
for |n — ml = 1 since under normal conditions the relation 
<C 1 B is fulfilled. To evaluate b, we now apply the 
approximation 



1 



1 



1 - 



1 



Xn 



(3.4) 



Because of the condition (2.13) the above approximation 
i s ju stified for all terms of the sum in the nominator of Eq. 
(|3.3[), since in these terms |n — n| is of the order of magni- 



tude of V An^ or smaller, and since x/(l + X^) < ^/n for 
any val ue o f x- The second term in the square brackets 
of Eq. (3.4) therefore leads to a contribution to b that 
small in comparison to the time-independent contribu- 
tion of the first term. Replacing fT„_i_„(i) by its initial 
value cr„_i_„(0) = V^l^n n this small contribution and 
taking into account that X]^o '^^^ri,n = -(- An^, we 
obtain after minor algebra a time- independent ex press ion 
for b. The latter can be substituted into Eq. ( 2.1?: ) to 
yield the approximative expression 



Aw = 7 



rx 



2(1 + xn) 



X 



1 



xn 



An^ 



(3.5) 



for the laser linewidth. 

In the next step we eliminate the injection rate r 
by m aking use of the photon-number balance equation 
( ^.12 ). For the case of the laser the latter takes the form 



rx sr^ 



2 ^ 1 

n=0 



1 



x(" + 1) 



Pn 



"f{n - Ub). 



(3.6) 



Here again Eqs . (3.1) and (3^) have been used together 
with Eq. (2.16), and the index of summation has been 
changed appropriately. We proceed by applying the same 
approximation scheme that led to Eq. (3.5) and perform 
the expansion 



1 



x(i 



n) 



I + xn 



1 + xin + 1) I + xn 
Using the condition ( p. 13 ) and the relation 
m ^ n, 



we obtain from Eq. 
the approximation 



rx 



2(1 + xn) 



1 



(3.7) 



(3.8) 



&) after simple transformations 



l+Ub 



xn 



(3.9) 



Whe n on ly the first term in the square brackets is kept, 
Eq. (3^) yields the familiar relation 



27 



1 

X 



1 

X V7 



- 1 



(3.10) 



where we have introduced the linear gain a = rx/2. 
However, in order to study the quantum limit of the 
linewidth, also the terms that are small in comparison 
to the leading term have to be taken into account. When 
we substitute the expression (|3.9| ) for the factor in front of 
the square brackets in Eq. ( |3.5D , keeping only the terms 
of lowest order in 1/n and x/(l + Xn), respectively, the 
contributions containing the variance cancel and we fi- 
nally arrive at the laser linewidth 



Auj 



1_ 
2n 



1 



— -I- 2nb 



xn 



(3.11) 



The quantum origin of the intrinsic laser linewidth is re- 
vealed by noticing that along our lines we would have 
obtained the result Auj = if the terms of the order 1 /n 
arising from the commutation relation for t he fi eld opera- 
tors had been neglected in Eqs. (3^) and (3^). The two 
limiting cases xn 3> 1 and xn <^ 1 correspond to a laser 
that is operated far above threshold or n ear threshold, 
respectively, as can be seen from Eq. (3.10). Therefore it 
is apparent that the linewidth becomes inversely propor- 
tional to the mean photon number n only near threshold 
and far above threshold. In the intermediate regions the 

dependence on n is more complicated. 

We note that because of the relation ( ^.1C| ) the 
linewidth can be expressed in terms of any three of the 
four paramete rs tT, y, 7, and a (or r, respectively). By 
inserting Eq. ( [3.10| ) into Eq. ( 3.11 ) and thus eliminating 



n, we obtain another useful expression for the linewidth. 



IX _ 
2 a 



7 



7 
a 



2nb 



(3.12) 



Eq. (3.11) describes the dependence of the linewidt h of a 



given laser on the mean photon number, while Eq. (3.12) 
describes the dependence on the gain or, alternatively, on 
the above-threshold ratio, defined simply as the normal- 
ized gain, a/7. Since the linear gain is easily measurable, 
this latter equation is the most important of the possi- 
ble expressions for the intrinsic laser linewidth from the 
point of view of experimental accessibility. 

In order to facilitate later comparison with the stan- 
dard result delivered by the phase-diffusio n mo del, we 
eliminate x from the linewidth expression (3.11). Using 
Eq. (3.10) one more times, we arrive immediately at 



a + 77 7 

Auj = — — h -Ub. 

In a n 



(3.13) 



Although it might seem from this expression that Auj is 
proportional to in the entire region of laser operation, 
this is not true since the gain a is not constant for a given 
laser but depends on the pumping rate and is connected 
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with n via the relation ( 3.10| ). In contrast to this, the 
cavity decay constant 7 and the saturation parameter x 
are fixed, being the saturation photon number for 
the lasing transition between the atomic energy levels. 

It is interesting to consider two important limiting 
cases. In the far- above-thresho ld lim it, where xn ^ 1 
or a 3> 7, respectively, Eqs. (3.11) - (B.13) yield the 
limiting value 



A..u„. = ^(l 



2nh 



la 



2nh 



(3.14) 



On the other hand, in the vicinity of the threshold, where 
0<q;/7— 1 <C 1 jl^, the linewidth depends on rl in a 
different way, described by Awthr = =(1 + J^fa)- This 
difference can be interpreted to be due to the fact that 
the contribution of intensity fluctuations to the linewidth 
increases when the above-threshold ratio decreases, as we 
shall show next. 

In a quantized description of the radiation field, in- 
tensity fluctuations are revealed in an enhancement of 
the photon-number variance as compared to the Poisso- 
nian value of An^ = n that corresponds to a constant 
intensity. Therefore, we first calculate the steady-state 
photon- number variance An^ of the laser, taking into 
account the presence of thermal photons. To do so, we 
apply an approximation method that has been developed 
previously by one of the authors in order to investigate 
the photon statistics in saturated multi-photon atom- 
field interaction We start from the steady-state 
equation 



dt 



(An^) 



(n - 2mi)p„ = 0, 



(3.15) 



and use the right-ha nd s ide of the m aster equation (2.5), 
together with Eqs. ( ^.2[ ) and ( ^.16 | ), i n orde r to express 
p„ „. By taking into account Eqs. ([3.l|) and (3^) and by 
appropriately changing the index of summation in the in- 
dividual terms, the resulting equation can be transformed 
to yield 



(n + l)(2n+ 1 - 2n) 



1 + X('^+1) 



2n'y{nb — n) 



n=0 

+7K(n + l)(2n + 1) - (1 + nb)n{2n - 1)]} = 0, (3.16) 

with a ~ rx/2. The following treatment again relics on 
the assumption of a strongly peake d pho ton-number dis- 
tribution subject to the conditions ( 2.13 ) that imply that 
the approxi mation ( |3.7D is valid. When the latter is ap- 
plied in Eq. (3.16), we obtain a term that is proportional 
=o"^f'n,n — Applying the approximation 



I + X^ 



2An2 
l + Xn 



7[2An^-n(l + 2n6)]. (3.18) 



Here again the relation = n^-f An^ has been taken into 
account, and small contributions have been neglected. 
Finally, we again make use of the lowest-order photon- 
number balance equation a/{l + xn) = 7 (cf. Eq.(3.10)). 
Thus, from Eq. ( |3.18 ), after simple transformations, we 
arrive at the relative photon-number variance of the laser. 



An2 / 1 , 

n \ xn, 



a — 7 



-(1 + n,). (3.19) 



Clearly, when the thermal photon number rib is not small 
in comparison to 1, its influence on the photon- number 
variance of the laser is crucial even for Ub <C n, as is its 
influence on the intrinsic laser linewidth. For Ub = Eq. 



(3.19) corresponds to the standard result that is known 
from the literature We note that in general the 

relative strength of the intensity fluctuations is charac- 
terized by the normalized quantity 



(at a^) - (ata)2 _ n{n - I) _ ^ _ Q 



n2 



(3.20) 



Here we introduced the Mandel Q-parameter 

^_An2 l + Ub j + Uba 
Q = — 1 = nh H = — — , (3.21) 



xn 



a ~ '-f 



where Eq. ( 3.191 ) has been applied. Since the conditions 
( ^.13 ) imply that Q <C n, the intensity fluctuations char- 
acterized by Eq. ( 3.20 ) arc extremely small under these 
conditions, and they vanish in the limit n —i- 00. There- 
fore intensity fluctuations can be considered to be a true 
quantum effect in the above-threshold regime of the laser. 
Nevertheless their influence on the laser linewidth cannot 
be neglected, because the non-zero value of the latter it- 
self is a small quantum effect only. 

Wi th the h elp of Eq. (3.21) the linewidth equations 
( |3.1l| ) - ( 3.13 ) can be finally cast into the form 



Alo 



1_ 
2n 



Q 



i + Q 



(l + rifc). 



(3.22) 



Obviously, close to the threshold, where a/7 — 1 ^ 1 
jlTf and hence Q ^ 1, the contribution of the intensity 
fluctuations to the intrinsic laser linewidth is the same as 
that of the phase fluctuations. In this case the linewidth 
is twice as large as it would be without intensity fluctu- 
ations, i. e. for (5 = 0. 



IV. DISCUSSION 



[n + (n — n)]^ ^ n + 3nAn , 



(3.17) 



which is justified because of the condition (2.13), we ar- 
rive at the equation 



Before discussing the analytic results in more detail, it 
seems appropriate to say a few words about their range 
of validity, determined by the applicability of ou r ba- 
sic assumptions ( 2.13 ). Making use of Eq. ( 3.19| ), the 
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second of the inequalities ( 2.13| ) can be transformed to 
yield the condition [(1 + nb)a/{a — 7)]^/^ <C n^/^. When 
the laser is operated e. g. 10% above threshold, i. e. 
for a/7 ~ 1.1, this requirement can be assumed to be 
fulfilled if n > lO'^ which, because of Eq. (B.IO), corre- 
sponds to X < 10^^. At higher above-threshold ratios 
our approximation is valid for even smaller mean photon 
numbers, or larger values of x, respectively. In general, 
because of Eqs. ( p^ and (pilgj) we can combine the 
two inequalities of the condition ( 2.13 ) to yield the re- 
quirement 



Vx(l + "6) < - - 1 < -, 

7 X 



(4.1) 



which has to be fulfilled for the linewidth formulae (3.11) 
- ( 3.13 ) to be vahd. Obviously, if x ~ 10~^ and n;, <C 1 
as in typical continuous-wave gas lasers, our results are 
already approximately valid when the laser is operated 
only more than 1% above t hresh old, i. e. for a/7 > 1.01. 

The linewidth formulae ( 3.11 ) - ( 3.13 ) constitute the 
main result of this paper. They differ from the standard 
laser linewidth formula, 



pd 



1 + a 7X a + 7 



4n 



7 



J_ 

2n' 



(4.2) 



which has been derived for n;, = in the so-called 
phase-diffusion model neglecting intensity fluctuations 
Ip^ . Here, for the last two steps, we used Eq. ( 3.10| ) 
and the equivalent relation a = 7(1-1- x^) order to 
transform the standard result. 

In Fig. 1 the linewidth is plotted for different operat- 
ing regimes of the laser. We emphasize that the curves 
representing our result, ( 3.12 ), are in perfect agreement 
with the numerical results, found previously for the laser 
linewidth by computing the two-time correlation func- 
tion of t he fi eld p^ . From a comparison of Eq. ( |4.2| ) 
to Eq. ( 3.13 ), it is obvious that the phase-diffusion re- 
sult is only a good approximation for the linewidth when 
a/7 « 2, whereas it underestimates the linewidth closer 
to the threshold. Higher above threshold, o n th e other 
hand, the linewidth is overestimated by Eq. ( |4.2| ) which, 
for n ^ 00, yields the intensity- independent residual 
linewidth Aw = 7x/4, instead of Aoj — 0, to be expected 
in the classical limit. Moreover, we conclude from Fig. 1 
that the linewidth can be approximated by Awiim, given 
by Eq. ( 3.14 ), provided that a/7 ^ 5. The difference be- 
tween our result and the standard one in the far-above- 
threshold region shows that the phase-diffusion assump- 
tion of the standard treatment overestimates the contri- 
bution of phase fluctuations to the linewidth the more 
the higher above threshold the laser is operated, i. e. 
the stronger the effect of the nonlinearity stemming from 
the gain saturation. On the other hand, since intensity 
fluctuations are neglected, the standard treatment under- 
estimates the linewidth in the near threshold regime and 
corrections are necessary to incorporate the effect of the 



super-Poissonian photon statistics, as has been empha- 
sized already by Lu ||ll[]. As the above-threshold ratio 
increases and the intensity becomes more stabilized, the 
effect of this underestimation decreases. Fig. 1 suggests 
that the effects of overestimating phase fluctuations, on 
the one hand, and underestimating intensity fluctuations, 
on the other, just compensate approximately when the 
laser is operated around 100% above threshold. 




1.2 1.4 1.6 1.8 2 2.2 
above- threshold ratio 



2.4 




2 3 4 5 6 7 8 
above- threshold ratio 

FIG. 1. Normalized laser linewidth Ati;/(x7) versus the 
above-threshold ratio 0/7 for nt = 0. The full line cor- 
responds to our formula, Eq. ( ^.12| ). For comparison, the 
curves resulting from the sta nda rd expression Aojpd (dashed 
line, corresponding to Eq. (1.2)) and from t he ap proxima- 
tion AtJiini (dotted line, corresponding to Eq. ( 3.14[ )) are also 
displayed. 

In the following we shall discuss the reasons for the 
discrepancy between our result and the standard one 
in more detail. For this purpose, we first consider the 
different approximation methods that are employed in 
the literature for the derivation of Eq. (4.2). In the 
most common approach, the master equation of the den- 
sity operator is transformed into an equation for its P- 
representation. After changing to polar coordinates by 
writing the complex field amplitude as e = rexp(i(/)), the 
exact evolution equation for the quasi-probability den- 
sity P(r, (j)) contains derivatives with respect to r and 
(f) to all orders Q. This is due to the nonlinearity of 
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the underlying master e qua tion, w hich is revealed by the 
denominators in Eqs. and (3^). In the standard 



treatment it is assumed that P does not change along the 
radial coordinate, corresponding to a neglect of intensity 
fluctuations, and that only derivatives up to second order 
have to be taken into account. The latter assumption is 
necessary to ensure phase diffusion, but it becomes less 
and less justifiable as gain saturation, hence nonlinearity 
of the laser equations, increases. In this case the simple 
phase diffusion model can no longer be applied. In addi- 
tion, mixed terms also become important in the evolution 
equation for P, which are products of a differential op- 
erator acting on the amplitude r and another one acting 
on the phase. These, in turn, lead to cross correlations 
between the fluctuations of the intensity and the phase. 

In the frame of the photon-number representation of 
the density operator, the standard result has been de- 
rived by applying the quantum fluctuation-regression 
theorem and approximately investigating the decay of an 
initial value of the field instead of the two-time correla- 
tion function. This is done by means of determining the 
lowest eigenvalue that characterizes a single decay rate 
for all non-diagonal density-matrix elements Pn.n-i [§• 
In a more precise treatment a single decay rate would 
have to be determined for the quantities y/np„^n-i [ p^ , 
since the average field follows from performing the sum 
over these quantities. Moreover, with increasing above- 
threshold ratio the influence of the nonlinearity also in- 
creases and therefore the quantum fluctuation-regression 
theorem cannot be applied anymore, in general ||20| . 

With respect to the Heisenberg-Langevin approach, 
we mention that a nonlinear c-number Langevin equa- 
tion can be derived for the complex field amplitude 
e = r exp(i0). The coupling of the fluctuations of the real 
amplitude r and the phase (p is clearly obvious from this 
Langevin equation. The derivation of the phase diffusion 
result, Eq. rests on implicitly making a factor- 

ization assumption for expectation values containing the 
complex field amplitudes in the denominator and their 
noise operators in the nominator Q]. Because of the in- 
tensity fluctuations near threshold, and because of the 
nonlinearity due to gain saturation far above threshold, 
it would be extremely difficult to go beyond this approx- 
imation. 

We finally conclude that the laser linewidth cannot be 
explained satisfactorily with the help of the simple as- 
sumption that the intensity is constant and the electric 
field phasor executes a random walk in the complex plane 
as described by phase diffusion. In the linear approxi- 
mation, valid near threshold, it is true that the behav- 
ior of the phase ffuctuations alone can be described by 
phase diffusion, but intensity fluctuations contribute to 
the linewidth as well. Farther above threshold, on the 
other hand, the simple model breaks down because the 
ffuctuations of the intensity and of the phase are coupled 
due to the nonlinearity of the gain, and the behavior of 
the phase fluctuations cannot be characterized as sim- 
ple phase diffusion. The frequently used procedure of 



disregarding the microscopic intensity fluctuations and 
assuming that the photon-number distribution is strictly 
Poissonian, on the one hand, and not properly taking 
into account or even completely neglecting gain satura- 
tion, on the other, in general does not yield sufficiently 
accurate results for the quantum-limited linewidth. 



V. CONCLUSIONS 

In this paper we have studied the quantum-limited 
linewidth of a good-cavity laser by determining the first- 
order correlation function of the laser field at steady 
state. It is the decay of this correlation function and 
not the phase fiuctuations alone that determines coher- 
ence properties such as, e. g., the visibility of interference 
fringes. By taking the Fourier transform we obtained the 
power spectrum of the laser field and derived an ana- 
lytical expression for the quantum-limited linewidth as a 
function of the mean photon number (see Eq. (3.11 )) o r 
of the above-threshold ratio of the laser (see Eq. ( 3.12 ). 
Our analytical result is in perfect agreement with the 
results of earlier numerical studies |1J]. We explicitly 
demonstrated the effect of a super-Poissonian photon 
statistics (see Eq. ( ^.22| )) and showed that near thresh- 
old the linewidth is considerably larger than the standard 
phase-diffusion result, (cf. Eq. (|4.2|)), where the inten- 
sity is assumed to be constant, or the photon statistics to 
be Poissonian, respectively. Although in most practical 
cases the laser linewidth is limited by the much larger 
technical noise and the intrinsic quantum limit cannot 
be reached, there exists a variety of proposals to reach 
or even go beyond the quantum limit with the help of 
sophisticated methods [^,0. Our results show that for 
a precise quantum mechanical description of the laser 
linewidth it is necessary to directly calculate the first- 
order correlation function of the laser field, thus implic- 
itly incorporating intensity fiuctuations as well as phase 
fiuctuations, and to properly take into account the non- 
linearity of the gain. 
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